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The neutron induced reactions on stable Carbon, Nitrogen, and Oxygen isotopes are investigated by using 
the Talys1.4 toolkit with the default parameters. The neutron incident energy covers a range from 0.20 MeV to 
85.00 MeV. For "C and '4N, the Talys1.4 results agree with the experimental data, while the parameters should 
be adjusted for 60. Some E, windows are found by comparing the main channels of n + C/N/O reactions, 
which induce element change. In these En windows, a specific element is activated to a different one while 
leaving the other element atoms unchanged. The results will facilitate the research of doping effects in organic 


materials by using neutron activation technique. 
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I. INTRODUCTION 


Doping is an important method to change the properties of 
materials. For example, the substitutions of some atoms in or- 
ganic semi-conductors promote the charge transport and sta- 
bility greatly [1, 2]. The general method for doping is chemi- 
cal synthesis which introduces new elements into the original 
material. Actually, neutron activation is an important method 
to make doping in material, which can even change one ele- 
ment to other one, directly or after the decay of unstable iso- 
topes produced by activation. In this situation, different to the 
chemical doping, different doping effects can be expected. 

Itis well known that for different nuclei, there are different 
threshold energy values for specific channels in neutron in- 
duced reactions [3]. For material having different elements, 
the different threshold energy values for neutron reactions 
provide a chance for specific activation if the neutron energy 
is selected in experiments. 

In this article, the neutron induced reactions on stable car- 
bon, nitrogen, and oxygen isotopes will be investigated by us- 
ing a Talys toolkit. The energy windows between the different 
isotopes will be analyzed. The theory is briefly discussed in 
Sec. II. The results are discussed in Sec. III, and a summary 
is presented in Sec. IV. 


Il. METHODS 


The optical model can describe the neutron induced reac- 
tion well in a wide range of incident energies. In Talys1.4, the 
ECIS-06 is implanted as a subroutine to deal with the optical 
model calculations [4]. The description of Talys1.4 and the 
implanted functions can be found in the user manual [5, 6]. 
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In statistical models for predicting cross sections, nuclear 
level densities are used at excitation energies where discrete 
level information is not available or incomplete. Several mod- 
els are implanted to describe the level density in Talys, which 
range from phenomenological analytical expressions to tab- 
ulated level densities derived from microscopic models. The 
constant temperature and Fermi-gas model is set as the de- 
fault parameter of level density at low excitation energy re- 
gion, while the Fermi-gas model is used in the high excitation 
energy region [6]. 

Due to the threshold energy exists for different channels in 
the neutron induced reaction on a nucleus, one channel can 
only happen if the incident energy of the neutron is above the 
threshold, which potentially provides energy windows only 
one channel can happen while the other channels are prohib- 
ited. In this work, the incident energy of the calculated re- 
actions ranges from 0.20 MeV to 85.00 MeV, and the default 
parameters in Talys1.4 are adopted. The main reaction chan- 
nels include (n, np), (n, p), (n, œ), (n, 2n) and (n, y). The 
calculated results are compared to the measured data, which 
are extracted from the EXFOR library provided by the Na- 
tional Nuclear Data Center (NNDC) [7]. All the natural abun- 
dance data and the half-life time of isotopes are taken from 
Wikipedia [8]. 


IIl. RESULTS AND DISCUSSION 
A. Neutron induced reactions on carbon isotopes 


The C element has two stable isotopes, '?C and PC, with 
natural abundances of 98.93% and 1.07%, respectively. Only 
the n + 12:13C reactions are calculated. 


1. n+ PC reactions 


The ?C(n, p)?B, "Cí(n, o)?Be, and "C(n, 2n)!!C chan- 
nels have been measured previously. !?B is an unstable nu- 
cleus, which decays to 12C via electron emission with a half- 
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life time of 20.20 ms. ?Be is a stable isotope. !!C is also un- 
stable, which decays to !!B via positron emission with a rela- 
tive long half-life time of 20.334 min. Thus the "C(n, 2n)! C 
and "Cín, np)!!B channels are the main channels which re- 
sult in element changes with the same final products, !!B. 
The results are plotted in Fig. 1. 
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Fig. 1. (Color online) The yield of residues produced in n + '?C 
reactions. The measured results are denoted by symbols, and the 
calculated results are represented by lines. 


First, for the Cn, o)?Be channel, the Talys1.4 re- 
sults are consistent with the measured results when E, < 
11.00 MeV [9], while the Talys1.4 results differ largely with 
the results measured by Stevens et al. [10] in the range of 
18.606 MeV < E, < 21.46MeV. In the energy range of 
20.00MeV < E, < 60.70MeV, the Talys1.4 results are 
also much larger than the calculated results by Dimbylow et 
al. [11]. Above E, > 11.00 MeV, the results of the "C(n, 
o)?Be channel have large difference, which suggests that fur- 
ther experiment should be performed for a systematic under- 
standing. 

Second, for the '*C(n, p) B channel, the measured re- 
sults by Kreger et al. [12] and Rimmer et al. [13] coin- 
cide when E, « 16.00 MeV, but differ when Ej is higher. 
The Talys1.4 results agree well with the measured results 
when E, « 16.00MeV, but overestimate the measured re- 
sults by Kreger et al. [12] in the range of 16.00MeV < 
E, < 22.00MeV. The calculated results of '*C(n, p)!?B 
in the range of 20.00MeV < E, < 60.70MeV by Dim- 
bylow et al. [11], which also uses the optical model, prefer 
the measured results by Kreger et al. [12]. Meanwhile, the 
Talys1.4 predicts similar results between "C(n, p)?B and 
Cn, np)!!B reactions when E, < 20.00 MeV. 

Third, for the !2C(n, 2n)!!C channel, which have been mea- 
sured by many groups, the measured data agree well when 
Fi, < 27.00 MeV. The results can be divided into two groups 
in the range of 27.00MeV < E, < 40.00MeV, of which 
the upper group was measured by Welch et al. [14], An- 
ders et al. [15], and Kim et al. [16]; and the bottom group 
by Brill et al. [17], Uno et al. [18], Brolley et al. [19], and 
Soewarsono et al. [20]. The calculated results by Dimby- 
low et al. [11] prefer the upper group results when En < 
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35.00 MeV. Kim et al. [16] measured the results in the en- 
ergy range from 55.00 MeV to 64.00 MeV. The calculated 
results by Talys1.4 largely overestimate the measured ones 
when En < 50.00 MeV, but agree with the measured results 
by Kim et al. [16]. The calculated results by Dimbylow et 
al. [11] underestimate the measured results by Kim et al. [16]. 

The threshold energies of the ?C(n, a)°Be, "C(n, p)" B, 
Cm, np)!!B and ?C(n, 2n)!!C channels increase, which are 
about 6.18 MeV, 13.64 MeV, 14.89 MeV and 20.30 MeV, re- 
spectively. The "C(n, y)'?C channel happens in the whole 
energy range, but have a much smaller probability compared 
to other channels. 


2. n4 PC reactions 
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Fig. 2. (Color online) The yield of residues produced in the different 
channels of the n + '?C reactions calculated by Talys1.4. 


No measured data for the neutron induced reactions on 
IC is found. The Talys1.4 calculated results are plotted in 
Fig. 2. The threshold energy values increase with the (n, œ), 
(n, 2n), (n, p) and (n, np) channels, with the values 4.13 MeV, 
5.33 MeV, 13.64 MeV and 16.50 MeV, respectively. The (n, 
y) channel has no lowest energy threshold. The (n, y) and 
(n, p) channels have relatively low probabilities compared 
to the other channels. In the whole energy range calculated, 
when E, < 15.00 MeV, the main channel is PC(n, «)!°Be, 
in which !°Be decays to !°B by electron emission with a very 
long half-life time 1.39 x 10° years; when E, > 15.00 MeV, 
the dominant channel is ?C(n, np)? B, in which ?B mainly 
decays to C with a half-life time 20.20 ms. 


B. Neutron induced reactions on nitrogen isotopes 


The N element has two stable isotopes, '4N and ?N, with a 
natural abundance of 99.64% and 0.36%, respectively. Only 
the n + !4:15N reactions are calculated. 
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Fig. 3. (Color online) The yield of residues produced in the different channel of n 4- 14N reactions: (a) "N(n, o)!!B; (b) '^N(n, pc; (c) 
"N(n, 2n) PN; and (d) NG, np) PC and “N(n, yy N. The calculated results are plotted as different lines. 


1. n4 "N reactions 


The calculated channels for the n+ !4N reaction are '4N(n, 
np) ?C, '^N(n, p)'^C, '^N(n, o)! B, ^N(n, y) PN, and N(n, 
2n)N reactions, respectively. The nuclei PC, PN, and !!B 
are stable, while ^C and '?N are unstable. !?N decays to ^C 
by positron emission with a half-life time 9.965 min. The (n, 
np), (n, p), (n, 2n), and (n, «) channels are the main ones 
which will result in element changes. 


The results of the n+ '4N reaction are plotted in Fig. 3. For 
clarity, the results are plotted in different panels. The '4N(n, 
o)! B and '4N(n, p)'^C have small threshold energies, which 
are very similar. The !^N(n, «)!!B channel has been mea- 
sured by different groups [21-25], and the '^N(n, p)'^C has 
also been measured [21, 23—25]. For the two channels, the 
data measured by different groups is consistent. The Talys1.4 
results agree with the measured data in the low incident en- 
ergies, but overestimate the measured results when the Æp in- 
creases. 


For the '^N(n, 2n) ^N channel, the measured results [26— 
30] are consistent when E, < 19.00MeV. When E, > 
24.00 MeV, the measured results by Brill et al. [17] and 
Yashima et al. [31] are relatively consistent, and the calcu- 
lated results by Dimbylow et al. [11] also agree with the 
measured data but have relatively large errors. The predicted 
threshold energy value by Talys1.4 is E, œ~ 11.00 MeV, but 
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Fig. 4. (Color online) The yield of residues produced in the different 
channel of the n + '?N reactions calculated by Talys1.4. 


the Talys1.4 cross sections overestimate the measured results 
in the whole energy range, which increases fast with E, and 
reaches maximum at E, = 24.00 MeV and decreases with 
E, when E, > 24.00MeV. Since PN decays to 13C, the 
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Fig. 5. (Color online) The yield of residues produced in different channels of n + !9O reaction. (a) O(n, «)?C; (b) '°O(n, p)!N; (c) O(n, 
Qn) O; and (d) On, np) PN and On, y)"O. The results calculated by Talys1.4 are plotted as lines. 


final production is the same as the '4N(n, np)?C channel. 
The 'N(n, np) "C channel has a low threshold energy value 
of about 6.00 MeV, and the cross section of the channel in- 
crease fast with En, which peaks at about £j = 14.00 MeV. 

The threshold energies increase in the order of (n, œ), 
(n, np), and (n, 2n), which are about 0.13 MeV, 5.67 MeV 
and 11.27 MeV, respectively. The (n, y) and (n, p) chan- 
nel happens in the whole £j range and forms a peak around 
19.00 MeVand 9.50 MeV with a wide width, but the cross 
sections are very small compared to other channels. 


2. n4 PN reactions 


No measured data for the n + !?N reaction is found. Only 
the Talys1.4 calculated results are plotted in Fig. 4. The cal- 
culated results for the (n, np) and (n, œ) channels have almost 
the same values when Ej < 12.00 MeV and the (n, np) chan- 
nel has much larger values than that of the (n, x) channel. At 
the same time, the cross sections of (n, p) and (n, 2n) chan- 
nels only have small difference when Ej > 17.00 MeV. The 
(n, np), (n, œ), (n, p), and (n, 2n) channels have relatively 
similar threshold energy values, which are about 8.57 MeV, 
8.09 MeV, 9.55 MeV and 11.56 MeV, respectively. The 
cross sections of the (n, y) channel are much smaller com- 


pared to the other channels. When En > 40.00 MeV, large 
fluctuations are found in the results of the (n, p) and (n, y) 
channels. 


C. Neutron induced reactions on oxygen isotopes 


The O element has three stable isotopes, '°O, 17O, and !80, 
with a natural abundance of 99.75%, 0.003896, and 0.205%. 
Only the n + 16:150) reactions are calculated. 


1. n+ 0O reactions 


In Fig. 5, the results of the n + !9O reactions are plot- 
ted. The channels include !6O(n, o) PC, !6O(n, p)!6N, !O(n, 
2n)O, !6O(n, np) ?N, and !6O (n, y) 1O. !6N decays to 
160 by emitting electron with a half-life time of 7.13s. Thus 
the main channels which make the element change are '°O(n, 
o) ?C and !6O(n, np) PN. 

For the !éO(n, o)?C channel, the Talys1.4 calculated 
results are consistent with these measured by Johnson ef 
al. [32], and Seitz et al. [33] except those by Divatia et al. [34] 
when Ep is smaller than 5.00 MeV. When En > 7.00 MeV, 
the measured results by Dickens et al. [35], Bormann et 
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Fig. 6. (Color online) The yields of residues produced in the 
different channels of the n + "O reactions calculated by Talys1.4. 


al. [36], and Dandy et al. [37] are consistent, but the Talys1.4 
calculated results are unable to reproduce measured data well. 

For the '°O(n, p)!ÓN channel, the measured results by Mar- 
tin et al. [38], Bormann et al. [39], and Seeman et al. [40] 
agree well. The measured results by Subashi et al. [41] and 
DeJuren et al. [42] also agree with those results but have rela- 
tively large difference. The Talys1.4 calculated results largely 
underestimate the measured results but have similar trend to 
the measured ones. 

For the '°O(n, 2n) ^O channel, the measured results by 
Yashima et al. [33] and Brill et al. [17] are in different en- 
ergy ranges, but for the overlapping range of En, the results 
have some difference. The Talys1.4 results overestimate the 
measured ones when E, < 40.00 MeV, but underestimate the 
measured ones when En > 40.00 MeV. 

e !6O(n, np) PN and !6O(n, y)!70 channels have not 
been measured. For the '°O(n, np)N channel, the proba- 
bility increases fast above the threshold energy of 10.57 MeV 
and has a maximum value around E, = 20.00MeV. The 
160(n, y) "O channel happens in the whole En range but have 
much smaller values. 

The Talys1.4 calculated threshold energies of the (n, o), 
(n, p), (n, np) and (n, 2n) channels for 160 are 2.36 MeV, 
10.25 MeV, 10.57 MeV and 16.65 MeV, with peaks form at 
around 10.00 MeV, 14.50 MeV, 26.00 MeV and 28.00 MeV, 
respectively. 


2. n4 "O reactions 


The main channels that the n + 180 reactions cover are 
180(n, a) 5c. 180(n, np)"N, 180(n, p) ^N, '80(n, 2n)!70, 
and '8O(n, y)!9O. UN can decay to !60 and 1O with a half- 
life time of 4.173 s; !5N can decay to !*O, C, or 17O with a 
half-life time of 622 ms, which are all stable nuclei (!^C has a 
very long half-life time). 
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In Fig. 6, the calculated results for the channels are 
plotted. The thresholds of the (n, œ), (n, 2n), (n, p), 
and (n, np) channels are 5.29 MeV, 8.49 MeV, 13.85 MeV 
and 14.49 MeV, and the peaks form at around 10.00 MeV, 
17.00 MeV, 18.00 MeV and 49.50 MeV, respectively. When 
Ea > 33.00 MeV, the yield of the (n, p) channel has large 
fluctuations with Fn. 


D. Comparison between main channels of n 4- C/N/O isotopes 


The differences between the threshold energies of the 
different channels make it possible to change one element 
to another. To change the element in organic materials 
specifically, the incident energy of the neutron should be se- 
lected to fit the energy window, as has been illustrated in 
the results above. The comparison between the thresholds 
of the main channel inducing element changes will show the 
En window more clearly. 

r '°C, both the final production of "C(n, np)!!B and 
12C(n, 2n)!!C are !!B since !!C decays to !!B via positron 
emission with a half-life time of 20.30 ms. The threshold 
energy of "C(n, «)?Be is about 6.18MeV. When E, > 
7.00 MeV, PC can be changed to ?Be. The threshold en- 
ergy of "C(n, np)!!B is about 14.89 MeV. When E, > 
16.00 MeV, C can be changed to ?Be and !!B. This pro- 
vides actual application of neutron activation on C isotopes. 
The chemically synthesized compound, in which a C atom 
is substituted by a B atom, demonstrates a novel molecu- 
lar engineering concept of organic semiconductors [43]. For 
^N, the most important channel is (n, 2n), which happens at 
very low neutron incident energy. Since the final yields in 
the n + 1^15N reactions are mainly carbon isotopes, we do 
not discuss the energy window for these channels. For !6O, 
the main channels are !°O(n, «)?C and !6O(n, np) N when 
E, < 15.00 MeV, with the production of °C and PN, respec- 
tively. There is also an £j, window between the two channels 
in the range from 4.00 MeV to 11.00 MeV. With larger En, 

e óO(n, 2n) ^O channel is opened. 15O decays to '4N by 
emitting a positron with a half-life time of 70.60 s, i.e., !O is 
changed to !^N finally. 

For a better understanding of the energy windows, the en- 
ergy above the Talys calculated threshold energy for each 
channel is plotted in Fig. 8. The numbers from 1 to 24 
represent the channels and the alphabet from a to w repre- 
sent the values of the threshold energies, of which are also 
listed as follows (the unit is MeV): 1, E, = O for the 
IN(n, p)'4C and (n, y) channel of C/N/O; 2, !^N(n, o)!!B 
(a = 0.13); 3, O(n, «)P?C (b = 2.36); 4, PC(n, o)!?Be 
(c = 4.13); 5, O(n, o)PC(d = 5.29); 6, ?C(n, 2n)?C 
(e = 5.33); 7, “NG, np)PC (f = 5.67); 8, "C(n,o)?Be 
(g = 6.18); 9, ^N(n, c)" B (h = 8.09); 10, !5O(n, 2n)"O 
(i = 8.49); 11, PN(n, np)^C (j = 8.57); 12, PN(n, p)PC 
(k = 9.55); 13, '°O(n, p)!6N (1 = 10.25); 14, !°O(n, np) PN 
(m = 10.57); 15, “Nn, 2n)?N (n = 11.27); 16, PN(n, 
2n) ^N (o = 11.56); 17, PC(n, p) PB (p = 13.635); 18, 
12C(n, p)PB (q = 13.645); 19, !5O(n, p)'5N (r = 13.85); 20, 
350(n, np) N (s = 14.49); 21, "C(n, np)!!B (t = 14.89); 22, 
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Fig. 7. (Color online) Comparison among the results of the main channels of n + C/N/O reactions which induce element changes. 
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Fig. 8. (Color online) The comparison among the threshold energies 
of different channels which induce element changes (See the text for 
explanation). 


3C(n, np)?B (u = 16.49); 23, O(n, 2n) ^O (v = 16.65); 
24, Cn, 2n)!!C (w = 20.30). Though the energy windows 
are clearly shown in Fig. 8, it should be carefully analyzed 
when the neutron is used to activate special compounds, in- 
cluding different elements. 


IV. CONCLUSION 


In this article, the neutron induced reactions on the stable 
C, N, and O isotopes are investigated by using the Talys1.4 
toolkit, which calculates the reactions in the framework of op- 
tical model. On the one hand, it is found that for !2C and !4N, 
the Talys1.4 results agree with the experimental data, while 
for '°O, the parameters in Talys1.4 should be adjusted for a 
better prediction. On the other hand, a systematic compari- 
son among the main channels of n + C/N/O reactions which 
induce element change are performed to find Ej windows 
among the original C, N, and O stable isotopes (by consider- 
ing the final production of the channel, i.e., direct change or 
indirect change through decay to a different element isotope). 
In these En windows, specific elements can be activated to a 
different one while leaving the other element unchanged. The 
results may help to study material modification by using neu- 
tron induced doping techniques such as in organic materials 
like the organic semiconductor. 
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